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1. Introduction 

Binary mixtures of Bose-Einstein condensates are of great interest due to their complex 
dynamical features and their role in the emergence of macroscopic quantum phenomena. 
Mixtures are usually made up of two species, that can also be hyperfme states of the 
same alkali atom pp. They generally display repulsive self-interaction and are confined 
by different potentials. Depending on the inter-species interaction, two classes of stable 
configurations are possible: mixed and separated. The latter are more interesting, since 
they allow the observation of phenomena such as symmetry breaking, e.g. in harmonic 
potentials, and macroscopic quantum tunnelling (2J [3]. Binary mixtures in harmonic 
traps have been investigated in a number of interesting experiments [U [5j |6] . 

Different approaches are possible in order to study the ground state of these systems. 
The binary mixture of two species of bosons can be rigorously described in a second- 
quantization formalism [7] . However, if the number of particles in the condensate is very 
large compared to the number of particles in the excited states, the fields associated to 
the two species can be treated as classical wave functions. This approach leads to the 
Gross-Pitaevskij equations |8J, that are obtained by minimizing the zero-temperature 
grand-canonical energy of the system. The ground state of the system can be thus 
determined by solving the Gross-Pitaevskij equation [5J QUI Ell EE2], or equivalently by 
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analytically or numerically minimizing the grand-canonical energy functional (3j EJ [13] . 
Analytical results are obtained only in particular cases, such as confinement by a hard 
wall trap |14| . harmonic or lattice potentials [15] and axisymmetric traps [IB] . A 
simplified approach is often used, based on the Thomas-Fermi (TF) approximation, 
that consists in neglecting the kinetic energy with respect to the self- and inter-species 
interaction energies [8]. This reduces the problem of finding the ground state of the 
binary mixture of condensates to a classical problem, related to the stability of a system 
of two interacting fluids. 

In this article we shall investigate this problem by adopting the following approach: 
given a system of two interacting condensates, confined in a generic external potential 
(that can be different for the two species), we will find general tools to determine the 
ground state of such system in the TF approximation. The article has the following 
structure. In Section [2] we set up the problem and introduce notation. In Section [3] 
we find a threshold value of the inter-species interaction parameter, above which mixed 
configurations cannot be the ground state of the system. In Section [4] we establish those 
conditions that determine which one of the possible separated configurations is the 
ground state. We conclude with an outlook in Section [5] Throughout this article, both 
cases of i) fixed numbers of particles and ii) fixed chemical potentials will be considered. 

2. Gross-Pitaevskij equations and Thomas-Fermi solutions 

We consider a system made up of two species of indistinguishable particles, labelled 1 and 
2, confined by generally different external potentials V\{x) and V^Oe). Self interaction 
and inter-species interaction are assumed to be repulsive. An example of such a system 
is a mixture of alkali atoms in two different hyperfine states [U E] • The two subsystems 
are described in a quantum field theoretical framework, by associating to each species 
the field operators ^i(x) and $2 0*0 • However, since we are searching for the zero- 
temperature ground state, we assume that all particles condense in the same wave 
function, and thus apply a Bogolubov shift [TF] and treat ipi{ x ) an d i>2{x) as classical 
fields, normalized to the average numbers of particles N\ and N 2 . The grand-canonical 
energy functional is 




(1) 



where 





k = 1,2, 



(2) 
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and m is the dimension of the system. By requiring that the energy is stationary, one 
obtains the coupled Gross-Pitaevskij (GP) equations 

h 2 



A + V x {x) + U n \Mx)\ 2 + U 12 \Mx)\ 2 ) M*) = ViMx) > ( 3 ) 
A + V 2 (x) + U 22 \Mx)\ 2 + U 12 \Mx)\ 2 ) fa(x) = p 2 ^ 2 {x) . (4) 



2m\ 
h 2 

2m 2 j 
In the above equations JJ\\ and XJ 22 are the self-interaction parameters between atoms 
of the same species, while U\ 2 is associated to inter-species interaction. Each of these 
parameters is assumed to be positive, since we are considering repulsive interactions. 
The solutions of ([3])-((4]) depend on the value of the chemical potentials p\ and p 2 , which 
are Lagrange multipliers. If pi and p 2 are fixed, the average particle numbers are free 
to vary. If, on the other hand, the particle numbers Ni and N 2 are fixed, the chemical 
potentials are chosen in such a way that the wave functions satisfy the normalization 
constraints A4 = Nk, for k — 1, 2. 

For the sake of simplicity, our analysis will be focused on one-dimensional systems, 
with the main results generalizable to higher dimensions. Moreover, it will be assumed 
that the potentials be continuously different iable, Vk € C 1 (IR). This class of potentials 
schematizes very well those used in trapping cold atoms. 

The Thomas-Fermi (TF) approximation, which will be applied in the following, 
consists in neglecting the kinetic energy contribution T to the energy functional Q. 
This approximation is justified if the number of particles is sufficiently high, since the 
self-energetic parts in Eq. (fil) are respectively O(iVf) and O(iVf), while the kinetic 
energy is 0{N 1 I ) + 0(A^ 2 ' ) |8J. As a consequence of the TF approximation, the 
grand-canonical energy becomes dependent only on the densities pi(x) = \ipi(x)\ 2 and 
p2(x) = l"^^)! 2 , an d will be indicated in the following as £tf(pi>P2)- 

Notice that, without loss of generality, one can reduce the analysis to the particular 
case U\\ = U 22 = 1. Indeed, by the scaling 



Pk ->■ Pkj \/Ukk, N k -» N k / \/Ukki V k ->■ Vk^/Ukk, Pk ->■ Pk \/Uk~k (5) 
one gets 

£tv(Px, Pi) = Ufa, p 2 ) - Pi (Af(pi) - Ni) - p 2 (Af(p 2 ) - N 2 ) , (6) 

with 

U (pi, P2) = ^ J (pl + p\ + 2otpip 2 ) dx + J (Vipt + V 2 p 2 ) dx, (7) 

and 

a = . (8) 

yUvJJ 22 

Incidentally, notice that the above reduction to a single parameter a applies also to the 
full energy functional Q, by scaling also the masses irik — > mk/y/Ukk- 

The critical points of the Thomas-Fermi grand-canonical energy functional are the 
solutions to the algebraic equations, 

px(x) + ap 2 (x) + Vi(x) = px, p 2 (x) + api(x) + V 2 (x) = p 2 (9) 
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and will be called the TF density profiles. Moreover, for fixed particle numbers N k , they 
are supplemented by the normalization conditions 

Pk dx = N k , (k = 1,2), (10) 

which fix the values of the chemical potentials p, k . 

In the following the supports of the TF densities p k will be denoted by S k . By 
assuming that a ^ 1, in Si 2 = Si H S 2 , where both condensates are present, the TF 
density profiles are 

I \ _ ^ ~ Vl ( X ) - « (A*2 - V 2 (x)) , . _ fX 2 ~ V 2 (x) -«(//!- Vi(x)) 

pi\ x ) — ; — 2 ' P2{ x ) — ; — o • y > 

1 — or 1 — cr 

In the regions S\\ = Si — S 2 and S 22 = S 2 — Si, occupied by only one of the two 
species, the solutions are respectively 

Pi(x)= f i 1 -V 1 (x), P2(x) = 0, (12) 

and 

p 2 (x) = p 2 -V 2 (x), pi(x) = 0. (13) 

The TF density profiles ( 1 1 )-( 13 ) are defined independently of the dimensionality of the 
system. 

Notice that the TF equations (|9| uniquely determine the functional dependence of 
the densities at a point x on the external potentials at the same point, the chemical 
potentials and the interaction parameters, once the supports Si and S 2 are given. On 
the other hand large freedom is left in the choice of the supports of the density profiles, 
for which uniqueness fails. Thus, extremely irregular configurations can be solutions of 
the TF equations. Among all possible solutions, one should pick up the minimizers. 

The rest of this paper will be devoted to deriving general rules for finding the 
minimizing configuration of the supports, in order to determine the ground state of the 
system, both if the numbers of particles or the chemical potentials are fixed. 



3. Mixed vs separated configurations 

The configurations of the binary mixture can be divided in two fundamental parts: 
separated and mixed. The TF densities are mixed in Si 2 — Si fl S 2 , where both species 
are present, and are separated in Su U S22 = £1 U S 2 — Si fl 5*2, where only one species is 
present at one time. A configuration is said to be separated if it does not contain mixed 
parts, and mixed otherwise. In this Section we will show that 

a th = 1 (U$ = y/UnU 22 ) (14) 

plays the role of a threshold value, above which separated configurations become 
energetically favored, both in the case of i) fixed numbers of particles and ii) fixed 
chemical potentials. This threshold holds independently of the particular external 
potentials V k (x). We proceed by treating separately cases i) and ii). 
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3.1. Solutions are confined 

We will first prove that under the assumption that the C 1 potentials are confining, that 
is 

Vk(x) — > +00, for \x\ —> 00, (15) 

with k = 1,2, all TF density profiles are compactly supported. We will see that this is 
a straight consequence of the positivity of the densities 

Pk{x) > 0. (16) 

We will prove that the supports Sk are bounded, by separately considering the sets 
with separated phases, Su = Si — S2 and S22 = S2 — Si, and that with mixed phases, 



S12 = SiD S 2 . By requiring that the solutions (12) and (13) be nonnegative we get 



S u C {x e M. I Vi(x) < = Vi \-oo, ml 

s 22 c{xeR\ v 2 (x) < m} = v 2 \-oo,ml (17) 



which are bounded by hypothesis. On the other hand, from ( 11 ) we get that every point 
x G Si 2 satisfies the conditions 

m - Vi(x) ~ a {m - V2(x)) > Q m - v 2 {x) -a(m- Vi(x)) > Q 

1 — a 2 ~ 1 — a 2 

They are easily proved to be equivalent to 

min{o;, a -1 } (m ~ V 2 (x)) < (m - Vi(x)) < max{a, a -1 } (m ~ V 2 (x)) , (19) 

which in turn imply that 

S12 c v^i-ocm] n v 2 -\-oo,m}, (20) 

so that S12 is compact. As a consequence Si = SnUSi 2 and S 2 = S22US12 are compact. 
3.2. Fixed numbers of particles 

If the numbers of particles Ni and N 2 are kept fixed, the chemical potentials are 
functionally dependent on the density profiles, since they have to be tuned in order 



to preserve the normalization conditions (10). The search for the zero-temperature 



ground state of the system reduces to the minimization of the TF grand-canonical 
energy functional that evaluated at the TF solutions reduces to the internal energy 
functional (|7j). 

3.2.1. Square well. A simple lemma will be now introduced (see e.g. [H]). Consider a 
binary mixture confined in an infinite square well, corresponding to a bounded interval 
S = [a, b] (with b > a) of the real axis. Let |S| be the (finite) length of the well. Since 
in this case Vk(x) = and the TF density profiles are flat, pk{x) = A^./|S|, the internal 
energy of the completely mixed configuration in S is 

Um= 2[S\ ^ + N * + 2aNlN2 ) ■ ( 21 ) 
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On the other hand, a separated configuration with Ni particles of the first species in a 
subset 5*i C S and N 2 particles of the second in S 2 = S— Si has densities pk{x) = N k /\Sk\ 
and internal energy 



IN 2 1 N 2 

"•( |S '» = 2^ + 2^ <22) 
which is in fact a function of the length |Si|. The minimum of U s is attained for supports 
Si and S 2 = S — Si such that 



Ni iV; 



2 



I'S'll \S 2 



(23) 



Condition (23) can be also expressed in terms of the densities in the separated 
configuration: 

Pi = 92- (24) 
Since the value of the internal energy at (|23|) is 



Us= 2\S\ ^ + N " + 2NlN2) ' (25) 
the minimizing separated configurations are energetically favorite if a > 1, while the 
mixed configurations are less energetic than all separated configurations if a < 1. For a 



binary mixture in a square well, this proves the role of Eq. ( 14 ) as a threshold value 



3.2.2. Selection principle and regularity of solutions. Notice that the minimizers (24) 



have a very high degeneracy, since every sets Si whose measures satisfy (23) correspond 
to possible TF configurations. Among them, despite the regularity of the potentials, 
there are extremely irregular configurations with highly entangled supports and infinitely 
many points of discontinuity {domain walls). However, such a phenomenon is a 
consequence of the TF approximation that, by neglecting the kinetic energy part T 
in ([I]), is also deprived of its regularizing effect on the densities. The kinetic term 
favors smooth density profiles. Indeed, due to the presence of the kinetic energy, the 
grand-canonical energy functional is defined on functions with square integrable 
(distribution) derivatives, and TF solutions are approximations thereof. In particular, 
in the 1-dimensional situation, a domain wall is a discontinuous approximation of 
a (absolutely) continuous function that changes between two values in a very short 
transition region with a large derivative. 

Thus, each domain wall of the TF solution would correspond to an additional cost, 
in terms of kinetic energy of the true solution, and the above-mentioned degeneracy 
would be lifted: among the TF degenerate minimizers, T would select those one(s) with 
the minimum number of domain walls. In the following we will make use of this selection 
principle and, in particular, we will only consider densities in the class of piecewise 
differentiable functions, pk G C 1 . That means that there is a finite subdivision of S& 
such that the restriction of pk to each subinterval [tj,tj + i] is continuously differentiable. 
Incidentally, notice that the potentials themselves can be assumed to be piecewise 
differentiable, without modifying our results. 
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Going back to the square- well case, by the selection principle, for a > 1, we 
end up with only two degenerate minimizers: one with S\ = [a, c], where c = 
(aN 2 + bN 1 )/(N 1 +N 2 ), and the other with St = [d,b], where d = (aN 1 + bN 2 )/(N 1 + N 2 ). 
Both configurations are separated and have a single domain wall. 



3.2.3. Generic potential. We now extend the above result to the case in which the 
mixture is not confined in a square well, but rather by generic continuously differentiable 
confining potentials Vk(x) with k — 1, 2. We will prove that, if a > 1, the ground state 
of the system cannot be a mixed configuration, and thus we can restrict our attention 
to the separated ones. 

Let the TF densities have a mixed configuration in Si 2 = SiC\S 2 , given by Eq. (11 ). 
Since the densities pk are assumed to be piecewise continuously differentiable, and their 
supports are compact, S 12 is the union of a finite number of compact intervals in which 
the pk are C 1 . Choose a segment u = [x , X\] of length \u\ = e > included in some of 
those intervals. The interval u contains 



n ; 



p k (x) dx = e(p k ), 



(A; = 1,2) 



(26) 



particles, where (•) denotes the average on oj. Since the potentials are C 1 , we can express 
them in u as 

V k (x) = V k (x ) + V^ k (x))(x - x ), (A; = 1,2) (27) 
for some G to. 

Since the pk are continuous in u, two points Xk exist in this segment, in which the 
functions equal their averages: 

Pk(xk) = (pk)- (28) 
Moreover, since the density functions are continuously differentiable in [xo, sci], taking 



into account (28) we can express them in each point of the segment as 
Pkix) = (p k ) + p' k (Vk(x))(x - x k ), 



(29) 



for some r]k(x) G u. The first derivatives appearing in (27)-(29) are all bounded in u. 
Thus, the internal energy of particles in the set u reads 



\(pi? + \(P2? + a( Pl )(p 2 ) + V 1 {x )(p 1 ) + V 2 {x )(p 2 ) 



We now divide u in two subintervals U\ = [xo, y] and u 2 = [y, xi] 
w = Wi U uj 2 . 



(30) 



(31) 



and replace the TF mixed densities in uj with flat and separated density profiles, 
preserving the particle numbers 

m 



Pi 
Pi 



\vi | 
n 2 

\u 2 \ 



for x G LUi, 
for x G u>2- 



(32) 
(33) 
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As a rule in choosing the bipartition of ou, we assume that the stationarity condition 
(23) for the internal energy in an infinite potential well is satisfied 

|w 2 1 _ n 2 

\iJl\ 7l\ 



(34) 



Taking into account the result (25), concerning the self-interaction and inter-species 
interaction parts, the potential energy of the set u with separated densities (32)-(33) 
can be expressed, after a straightforward manipulation, as 



Ms 



;(pi) 2 + 2<P2> 2 + (pi)(p 2 ) + V 1 (x )(pi) + V 2 (x )(p 2 ) 



(35) 



The net change in the total potential energy, due to the replacement of the mixed 
densities in [xo,xi] with the separated ones, is 

e(l-a)( Pl }( P2 ) + 0(E 2 ). (36) 



Us - u r 



For sufficiently small e, the sign of 8U is determined by the first term in (36), unless 
a — 1. 

If a > 1, the result SU < implies that, given a point of a mixed configuration, 
there always exists a neighborhood in which one can construct a separated configuration 
with lower energy. Since Si 2 is compact, we can find a finite subdivision Si 2 = 
[xo,xx] U . .. U [x n -i,x n ], such that the above construction can be performed in each 
segment Xj]. Thus, the minimizers are separated configurations, if the particle 

numbers are fixed. Analogously, one can show that if a < 1 the internal energy of a 
separated configuration is always larger than the energy of a mixed configuration with 
the same particle numbers. Thus, even in the case of varying external potentials the 
value (14) acts as a threshold between mixed and separated ground states. 



3. 3. Fixed chemical potentials 

If the chemical potentials \L\ and [i 2 are fixed, the average particle numbers are free to 
vary. In order to find the ground state of the system, one has to find the minimizers 
of the grand-canonical energy £ TF (pi, p 2 ). As in the case of fixed particle numbers, the 
elementary case of a binary mixture in an infinite square well will be analyzed first. 
Then, we will try and find general results in the case of different piecewise continuously 
differentiable confining potentials. 



3.3.1. Square well. It is clear from (11 )— ( 13 ) that fixing the chemical potentials 
corresponds to fixing the density functions, with the choice of the supports leading to 
different numbers of particles. We will neglect an inessential annoying constant in ^ 
by setting N x = N 2 = 0. 

In the simple case of an infinite square well S with V\ — V 2 — inside the well, the 
only values of the chemical potentials that have physical meaning are the positive ones, 
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as it emerges from (17) and (20). Let us first consider a completely mixed configuration 
in the well with density profiles 

m _ Pi - afJ-2 m _ P2 - an\ , 

Pi ~ 1 2 ' Pi ~ 1 2 ' v 

1 — or 1 — cr 

Such a configuration has a physical meaning if both densities are non negative. Thus, 
if a < 1 the numerators in (37) must be nonnegative, while if a > 1 the numerators 
must be nonpositive. Densities in a mixed configuration are not defined at the threshold 
value a = 1. Conditions on the positivity of the densities set a bound on the values of 
a which are compatible with the chosen chemical potentials: 

a (a 1 , a u ), 

where 



(38) 



a 



mm 



Pi P2 
P2' Pi 



< 1< a u :-- 



max 



Pi P2 
P2' Pi 



The problem of nonphysical values of a does not arise in the case 

P2 = Pi, 



(39) 



(40) 



which will prove to be a very relevant physical situation. If (40) holds, the boundaries 
(39) coincide, and the densities are well defined for all a. Note that condition (40) 
exactly corresponds to the minimum condition of the internal energy in the separated 
phase, since, by taking into account ( 12 )-( 13 ), it implies 

p\ = p s 2 . (41) 



If the solutions (37) are plugged in the definition of the grand canonical energy, it is 
possible to express it in terms of interaction parameters and chemical potentials 

Pi + P2 _ 2flifi 2 0i 



\S\ 



2(a 2 -r 



(42) 



If instead separated solutions are considered, with the first condensate confined in a 
region Si and the second in S2 = S — Si, the grand-canonical energy is a function of 
the length \Si\ and reads 



£s(\Sl\) = -\Si\y ~\S - Si|y 



Since (43) is linear in the length \S 
-\S\ max 



Pi P_| 
2 ' 2 



it is clear that its minimum value is 



(43) 



(44) 



Hence, if pi < pi, the minimum of the grand-canonical energy for separated 
configurations corresponds to Si = S, while if p 2 > Pi, it corresponds to S 2 = S. 
In both cases, the minimizer is in fact a single condensate configuration. Only if (40) 
holds, separated configurations are allowed. Moreover, their energy is stationary with 
respect to changes in the partition of S. 

The energy in the mixed configuration and that in the minimizing separated 
configuration will be now compared in detail, as a varies from zero to infinity, in the 
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cases fi2 < A*i and = A*i, since if /i2 > A*i the physical situation is specular with 
respect to the first case. We remark that the energy of each separated configuration is 
always independent of a. 

In the case ^ < if the inter-species interaction is absent, a = 0, the mixed 
configuration is favorite, since its energy £^ reads 

£l = -\S\(^ + f) <£ s . (45) 

If a increases, the energy of the mixed configuration grows, until it reaches a local 
maximum at or, which marks the beginning of the nonphysical region. At this point the 
density of the second species vanishes, and thus the energies of the mixed and separated 
configurations are equal: they are in fact both single-condensate configurations, whose 
energy is 

Sl = -\S\£ = E.. (46) 

For a 1 < a < a u , the separated ones (not only the minimizing one) are the only 
configurations that have a physical meaning. If a = a u , the mixed configurations 
become physical again, but their energy is higher than that of the single-condensate 
configuration 

Sl = -\S\^>£ S . (47) 

The single-condensate configuration remains energetically favored for a — > oo, since its 
energy is constant, while the energy of the mixed configuration (42) vanishes as a -1 . 

We now consider the case //2 = It was observed that this is the only case in which 
real separated states minimize the grand-canonical energy S s , and thus it is possible for 
such configurations to be the ground state of the system. Moreover, we stressed that 
if condition (40) holds, there is no nonphysical region for the mixed configurations, as 
a varies from zero to infinity. By plugging condition (40) in (42), we find that the 
grand-canonical energy in this case is never singular, and reads 

S m = -\S\-^. (48) 
1 + a 



If (48) is compared with (43), that for \i\ = fi 2 reads 

£ = -|£|Y (49) 

and is independent of the partition, mixed configurations are found to be favorite if 
a < 1, while separated configurations have smaller energy if a > 1. Thus, even in the 
case of fixed chemical potentials and infinite square well external potential, the value 
(14) proves to be the discriminant value between mixed and separated ground states. 
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3.3.2. Generic potential. Consider TF density profiles p k (x). In Su — Si — S2 we get 
(with the convention Ni = N 2 = 0) 

£ s (1) = \j fai - Vi? dx- [ (//a - Vr) 2 dx = -\ [ {p x - l^) 2 tfc, (50) 

and in S22 = S2 ~~ 



^ S " 2 



I / {^2-V 2 fdx. (51) 



Therefore 



with 



£* = —[ Vi(x) 2 dx-±f V 2 (x) 2 dx, (52) 

•/ Su J S22 



V k (x)=p k -V k (x), (k — 1,2). (53) 
In 5i2 = Si H S*2 we get 

£ m = - on 1 2 , / (v^x) 2 + V 2 (x) 2 - 2aV 1 (x)V 2 (x)) dx. (54) 
— a ) J Sl2 \ J 

Compare a mixed TF configuration in a set S with a configuration with only one species, 
say px, 

58 = £ s (1) -S ra = ] 2 jv 2 dx+ 2{l l _ a2) J (Vi + V 2 - 2aV 1 V 2 ) dx + 
' V 2 -aV^\ dx. (55) 



5 



2(1 - a 2 ) 

Thus, when a > 1 one gets 68 < and separated configurations are energetically 
favorite. Therefore, in this condition, we have to search for the ground state among the 
separated configurations, which is the aim of the next section. 

3-4- Remarks 

From the results obtained in this Section, it clearly emerges that, if a > 1, that is 
U12 > UI2, the ground state is in a separated configuration. Moreover, according to 



the selection principle introduced in Sec. 3.2.2, the following analysis will be restricted 



to piecewise continuously differentiable solutions of the TF equations, p k G C l , whose 
discontinuities can be due to the presence of a finite number of interfaces. 

4. Minimizing separated configurations 

In Section [3] it was shown that separated configurations are energetically favored if 

a > 1 , (56) 

for generic continuously differentiable confining potentials. We will now find a way to 
determine which of these configurations is the ground state of the system, and which 
can be regarded to be locally stable or unstable. We shall again discuss separately the 
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cases of fixed numbers of particles and fixed chemical potentials, underlining analogies 
and differences between them. 

When only piecewise continuously differentiable density profiles are considered, a 
separated configuration can be characterized by the property that the supports Si and 
S 2 do not intersect, except at a finite set of points. The intersection points correspond 
to a set of domain walls separating the first and the second species. Si and 5*2, being 
compact, are thus unions of intervals, which can be bounded by i) two domain walls, ii) 



a domain wall and a zero of the TF density profile, or iii) two zeros. From ( 12 )-( 13 ), 
the possible zeros Q k \ with j = 1, . . . , m and k = 1, 2, of the densities p k are subject to 
the condition 

m = V k {Cf ] ) . (57) 

In general, if we search the ground state among the separated configurations, we have 
to deal with the minimization of the sum of two decoupled functionals of the kind 

U(p 1 ,p 2 ) = V(p 1 ) + V(p 2 ). (58) 



We remark that, except for condition (56), that enables one to establish that the ground 
state is a separated configuration, the parameter a plays no role in the search for minimal 
separated configurations. 

We will proceed by fixing the number n of domain walls and determine the set 
of positions R = (R\, . . . ,R n ) for which the considered functional, internal or grand- 
canonical energy, is, at least locally, minimized. We will then compare the minima 
corresponding to different numbers of walls. 

4-1- Fixed numbers of particles 

In this case our aim is to find the stationary configurations of the internal energy 

U{ Pl ,p 2 ) = \ ! (pi(x) 2 + 2V l {x)p l {x))dx+ l - f (p 2 (x) 2 + 2V 2 (x)p 2 (x))dx(59) 

J Si J S2 

with respect to small variations of R, under the condition that the numbers of particles 
N k remain fixed, J Sk Pkdx = N^. Then, conditions will be set for these stationary 
configurations to be local minima. It can be easily inferred that the chemical potentials 



in ( 12 )-( 13 ), which are used as Lagrange multipliers to normalize the density profiles, 
and thus depend on the supports S\ and S 2 , are functions of the domain wall positions 
R. We re-express the separated TF density profiles, explicitly showing this additional 
dependence: 

p k (x;R)= p k (R)-V k (x), (k — 1,2). (60) 
The zeros C- > subject to condition (57), are also functions of R. Since the choice 



of the position of the domain walls completely defines the density profiles and their 
supports, the internal energy (59), evaluated at stationary TF densities, can be viewed 
as a function of R 

w(piM),p 2 (-;i2)) :=U(R). (61) 
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We now assign to each domain wall a dichotomic variable Sj\ Sj = +1 if it is the upper 
border of a interval containing the first species (and thus the lower border of an interval 
containing the second one), and —1 in the complementary case. 

By taking the first derivative of (61) with respect to a generic Rj and using 
normalization conditions, one gets 

dU(R) _ sj 



dRj 2 



(p 2 (R j ;R) 2 -p 1 (R j ;R) 2 ^) . (62) 



The stationarity condition is obtained by setting to zero the derivatives (62 ) with respect 
to all the positions of the n domain walls, yielding 

p 1 (R j ;R) = p 2 (R j ;R) Vj = l,...n. (63) 



These are clearly analogous to the minimum conditions (24) in the case of an infinite 
potential well. If the densities in (63 ) are expressed as functions of the external potentials 
and the chemical potentials, it becomes clear that the position of the domain walls in a 
stationary configuration are characterized by the fact that the potential 

<p(x) := V x {x) - V 2 {x) (64) 

is equal for all Rf 

<p(Rj) = m(R) - 1*2(3). (65) 

Let us consider the equation 

<p(x) = f, (66) 

with / a constant. The number of its solutions fixes the maximal number of walls in a 
stationary configuration. 

It is particularly interesting the case in which the external potentials for the two 
species are proportional, 

V 2 (x) = pV x {x) =: /3V{x) (67) 

with > 0. This happens, e.g. when the two condensates feel the same potential before 
the scaling (J5J), and in such a situation 



In the case of proportional potentials (67), if the equation V(x) = v, with v a constant, 
has n solutions, there cannot exist stationary configurations with more than n domain 
walls. In this case, the domain walls are placed at positions characterized by the same 
potential, which must be equal to 

V{R j ) = — , Vj = l,...,n. (69) 

As a consequence, the densities of the same species must be equal at the edge of each 
domain wall. The values of the densities at the edge of all domain walls in the case of 
proportional potentials will be indicated as p = px(Rj) = p 2 (Rj). 
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Figure 1. (Color online). Examples of maximal stationary configurations in a 
symmetric double- well potential. The blue dotted lines represent the external potential, 
in arbitrary units. The density profile of the first species is in dark grey (red in the 
online version), while the density profile of the second one is in light grey (green in the 
online version). In (a) the minima of the potential are occupied by the first species, 
which is the less self-interacting one. In (b) the minima are occupied by the second 
species. Computation of the second derivatives shows that only configuration (a) is 
(at least locally) stable. 



Henceforth, we shall call maximal stationary configurations those ones in which a 
domain wall is placed in each of the real solutions of ( 69 ) , except for the case in which 
one of the solutions is a stationary point for V(x). Two different examples of such 
configurations are shown in Figure [1} In the following, it will be shown that for (5 w 1 
the ground state of the system is usually in a maximal configuration. 

In order to determine if the stationary configurations are in fact minima of (61), 
the Hessian matrix H at the stationary solution has to be computed. By deriving (62) 
once more with respect to Rj, we find that d 2 U/dR"j has two contributions: the first 
one is related to the dependence of the external potentials on the point Rj, the second 
one to the dependence of the chemical potentials on the position of the domain walls 
R. On the other hand, if (62) is derived with respect to Rk with j ^ k, only the second 
one of the above mentioned contributions survives. By taking into account conditions 
(63), the diagonal elements of the Hessian matrix in the stationary configurations reads 

Hjj = s jPl {R f R)(p'(Rj) +LL + ^j Pl (R j; Rf , (70) 

while the non diagonal elements are 

Hjh = sjSk f jjj + J Px{Rj\ R)fh{Rh) R) ■ (71) 

For large numbers of particles, so that \S^\ are sufficiently large, all terms depending 
on the inverse length can be neglected, and conditions for the Hessian matrix to be 
positive definite, and hence for the corresponding stationary configuration to be a local 
minimum, are easy to find: 

(p'(Rj) > V((Rj) > V^Rj) if Sj = +1 , 

^(^•)<0 & V((Rj) < Vi{Rj) if s 3 = -l. (72) 

For smaller numbers of particles the complete Hessian matrix has to be diagonalized (e.g. 
in a numerical way). A simpler and relevant case is that of equal external potentials, 
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in which we have already remarked that in a stationary configuration the values of the 
densities of each species must be equal at all the domain walls. We obtain 

H% = 6 jkaj + (-lY +i C , (73) 
where we have taken into account that SjSk = (— l) 3+k and defined 

aj = Sjp(l -P)V'(Rj) (74) 
as the intensive and purely diagonal part, and 

c= {w\ + m)? (75) 

as the lenght-dependent term, which vanishes for large numbers of particles. In this 
limit, when /3 < 1 the condition for a stationary configuration to be locally stable is 
that the potential V(x) be increasing at all the R/s which are the upper (right) border 
of an interval containing particles of the first species (sj = +1), and decreasing at all 
the R/s which are the lower (left) border of an interval of the same kind (sj = — 1). 
Intuitively, the less self-interacting condensate [since we supposed that < 1 in Eq. (68)] 
tends to occupy regions of the real axis in which the potential is lower, while the most 
self-interacting one is pushed into regions where the external potential is higher. For 
small numbers of particles, such configurations continue to be minima for (61), since if 
the potential is increasing when Sj — +1 and decreasing when Sj = —1, in Eq. (73) 
is the sum of two positive definite matrices, and hence it is positive definite. Moreover, 
it is possible that even a configuration in which 

aj < for some j G {1, . . . , n} (76) 

becomes stable, which is impossible in the thermodynamic limit. However, bounds on 
the stability of such configurations can be found if some necessary conditions for (73) 
to be positive definite are tested. First, since |-£^ p | < (Hf? + H l -j)/2 for all pairs of 
indices, then 

ai + aj>0, Vi,j = l,...,n, (77) 

which implies that only one of the {%}, say aj, can be nonpositive for a stable 
configuration, and moreover 

\dj\ < minaj. (78) 

i¥=3 

On the other hand, by applying on the other hand the necessary condition det H > 0, 
another constraint can be established, 

C 

\a-A < == r , (79) 

with the upper bound vanishing in the thermodynamic limit. 

There are two kinds of maximal stationary configurations: the first one is 
characterized by the fact that the external potential at each point of Si is smaller than 
the potential at each point of S2, while the second one is characterized by the opposite 
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V(x) =v 



Figure 2. (Color online). Non-maximal configurations in a double-well potential. 
Lines and shades have the same meaning as in Figure [I] In (a) the first species 
occupies a region around the potential barrier, where the potential is higher than the 
value v in the domain walls. In (b) the second species occupies one of the minima of 
the external potential, where it is lower than the value v in the domain walls. 



situation. The latter, however, is not stable, since the diagonal part of its Hessian matrix 
contains Oj < for alii = 1, . . . ,n. We can thus limit our attention to the first kind 
of profiles, which we call the maximal stable configurations (see Figure [T]). We will now 
prove that if a locally stable configuration is not maximal, there are conditions ensuring 
that it cannot be the ground state of the system. The proof is based on the fact that in 
a non-maximal and locally stable configuration one of the following situations emerges: 
either the potential at a point of Si is greater than the potential at the domain walls, 
or the potential at a point of S% is smaller than its value at the domain walls. Examples 
of non-maximal configurations are represented in Figure [2] We remark that both these 
situations can be present in the same configuration. We start by considering the first 
one. Let 

v = V(Rj) (80) 

be the potential at the domain walls. Assume that the potential at a point Xq G Si be 
such that 

V(x )=:V 1 >v. (81) 

We exclude the case V± = n\, implying a vanishing density at x . (This case will be 
eventually considered as a limit.) We now consider a subinterval u e of S2 of length 
e > 0, with a domain wall as one of its edges, and an interval ui^ of length rj > 0, which 
is a neighborhood of xo, and impose that the number of first-species particles in uj v be 
equal to n and that of second-species particles in u e be n /(3. Since the potential is 
regular and the densities are supposed to be regular between each pair of domain walls, 
one gets 

n = VPifa) + 0( V 2 ) = ef3p 2 (R 3 ) + 0(e 2 ) . (82) 

This equality implies a relation between the length of the considered intervals, depending 
on the ratio of the densities, which can be expressed in terms of the potentials by using 



condition (63): 



Pi{xq) Hi - Vi 
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The potential energy of the two selected intervals is given by the sum of the contributions 



1 n 2 

«r = vno + + 0(e 2 ) , (84) 



,( 2 ) 



1 n. 



2 



4 ) = V 1 n + -^ ) +O(e>). (85) 

We now replace the original density profiles with flat density profiles which preserve the 
numbers of particles. In particular, we fill u t , which initially belonged to S% , with the 
first-species condensate with a density 

Pi = —, (86) 
and uj n with the second-species particles with a density 

With these new density profiles, the internal energy of the intervals becomes the sum of 
the terms 

«<" = tmo + \^ + 0(« 2 ) , (88) 

The total variation of the internal energy induced by this change is 

SU = «f > + - «?> - «f> = f - l) - |) + 0( E 2 ). (90) 

Since /3 < 1, we find that if 77(e) > e, i.e., up to first order in e, 



Ml - % 0' 

and it is always possible, for sufficiently small e, to find a density profile that preserves 
the numbers of particles, whose energy is smaller than the energy of a non-maximal 
stable configuration. It can be observed that if the limit V\ — > \L\ is taken, condition 
(91 ) is certainly satisfied. This means that in the ground-state configuration the intervals 
of the support Si of the less self-interacting species cannot be bordered by a zero. 

If the case in which there exists a point in xq G S2 where the external potential is 
lower than its value in the domain walls, 

V(x ) =:V 2 <v, (92) 

we find, with the same procedure as in the previous case, that if the following inequality 
is satisfied 

^<e, (93) 

P2 - V 2 

the considered non-maximal configuration can never be the ground state of the system. 
Since usually (3 2 = U\\jU%i ~ 1 (see, e.g., the hyperfine states of 87 Rb, |2]), conditions 
(91)-(93) set a very stringent limitation on the possibility that a non-maximal stable 
configuration be the ground state of the binary mixture. 



Binary BEC in a confining potential 



18 



4-2. Fixed chemical potentials 

The results in the case of fixed chemical potentials is very similar to that of fixed numbers 
of particles, in the thermodynamical limit. The functional to be minimized by separated 
configuration is (we set Nx = N 2 = 0) 

£ (Pi, Pi) = J Qp? + Vipi - faPi^J dx + J Qp^ + V2P2 - P2P2^ dx . (94) 

An important difference with respect to the previous case is that, since the chemical 
potentials are fixed and not subject to normalization conditions, the TF density 
functions are completely independent of the positions of the domain walls. Thus, 
the functional (94) depends on R only through the domains of integrations, which are 
determined by the supports of the density profiles, and it can be seen again as a function 
of the domain wall positions: 



£{ Pl (-,R),p 2 (-,R)) :=E(R). (95) 

The stationarity conditions are exactly the same as in the case of fixed numbers of 
particles, since the first derivative with respect to a generic Rj reads 

dE(R) _ sj 



dRj 2 



[p 2 {R, ] R) 2 - Pl {R j] R) 2 ') . (96) 



However, the Hessian matrix in the stationary configurations is diagonal, as in (70). The 



absence of the non-diagonal part lies in the fact that the first derivative (96) depends 
on the position of the domain walls only through the external potentials. The stability 
criterions for a stationary density profile are the same as in the case of fixed numbers 
of particle, if the thermodynamical limit is considered: a configuration is stable if and 



only if conditions (72) are satisfied for all j. If the two species lie in the same external 
potential, these conditions reduce to SjV'(Rj) > 0. 

Even in the case of fixed chemical potentials, it is possible to show that if the 



potentials are proportional (67), there are limitations on the possibility that a non- 
maximal stable configuration be the ground state. Indeed, if at the domain walls we 
have V(Rj) = v, it can be shown that if there exist a point in Si, where the potential 
is V\ > v and satisfies 

(97) 

then the grand-canonical energy of the configuration is higher than the energy of another 
configuration corresponding to the same chemical potentials. We observe that the 



stationarity condition (63) (with the chemical potentials independent of R) implies that 



the chemical potential fa be greater than fa. Thus, the left hand side of (97) is an 
increasing function of V\. Moreover, the equality is saturated for V\ = v. We remark 
that if the limit Vi — > fa is considered, that is if the density profile of the first species 



condensate has a zero, condition (97) is certainly verified, and then such a configuration 
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cannot be the ground state of the system. On the other hand, if there exist a point in 
S 2 where the potential is V 2 < v, satisfying 

V 2 



Mi 



> 13 



(98) 



- V 2 

then there exist another configuration corresponding to the same chemical potentials, 
which has a lower grand-canonical energy. In this second case, the left hand side of (98) 
is a decreasing function of V 2 , and the equality is again satisfied by V 2 = v. 

Here, we sketch the proof in the first case, the second case being analogous. Let 
us suppose that at a given point x G Si the potential satisfies V(xo) = V\ > v , and 
consider two intervals u and x °f the same length e, u lying in 5*2 and bordered by a 
domain wall, and % lying in S\ and containing xq. The chemical potentials are fixed, 
thus the functional form of the density functions does not depend on the positions of 
the domain walls. We now replace in u the first species with the second one, and in 
X the second species with the first one, using again TF density functions. Taking into 
account condition (63), we find that the difference between the energy of the final and 
the initial configurations is 







/i2 ~ Vl 



0(6 2 ) • 



(99) 



For sufficiently small e, the final configuration is energetically favored with respect to 
the initial one, if condition (97) applies. We finally note that, unlike in (91)-(93), the 
conditions (97)-(98) are independent of the value of v. Even in this case, since usually 
f3 2 = Uu/U 22 ~ 1, conditions (97)-(98) preclude any non-maximal stable configuration 
from being the ground state. 



5. Conclusion 



We have studied the Thomas-Fermi equations, for a system of two Bose-Einstein 
condensates confined in generic potentials. We have emphasized the role of the limiting 
value U\2 = VUnU 22 in determining if the ground state of the system is a mixed 
configuration or a stationary one, by assuming that the external potentials be regular. 
We then determined a set of conditions to be satisfied by locally stable separated 
configurations. Then we looked for the ground state among the possible locally stable 
configurations, and found that those with a maximal numbers of domain walls are 
usually energetically favorite. The results presented in this article enable us to find the 
ground state of binary mixtures in multi-well potentials, given either the numbers of 
particles or the chemical potentials. 

It would be interesting to analyze the changes that a correction to the TF 
approximation, including the kinetic energies, would introduce in such a picture. If 
the numbers of particles are sufficiently high, the TF approximation is very accurate. 
Nonetheless, TF density profiles correspond to diverging kinetic energy, due in particular 
to the discontinuities at the domain walls. The kinetic parts intervenes by regularizing 
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the TF solutions, at the expense of an increase in the potential energy, especially in a 
neighborhood of a domain wall. This could lead to an inversion in the energetic diagram, 
in which configurations with few domain walls could become energetically favorite with 
respect to maximal stable configurations. This inversion has already been numerically 
studied in the simple case of a harmonic potential [2J, but the tools introduced in this 
article uncover the possibility of extending this kind of analysis to generic multi-well 
potentials, such as arrays of optical traps, which are now within experimental reach. 
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